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ABSTRACT 

We consider the possibihty of constraining decaying dark matter by looking out through the Milky 
Way halo. Specifically we use Chandra blank sky observations to constrain the parameter space of 
sterile neutrinos. We find that a broad band in parameter space is still open, leaving the sterile 
neutrino as an excellent dark matter candidate. 

Subject headings: dark matter — elementary particles — neutrinos — X-rays: diffuse background 



1. INTRODUCTION 

The particle nature of the dark matter remains a mys- 
tery. This is contrasted with the firmly established dark 
matter abundance, which has been measured using the 
cosmic micro wave background and large scale struc ture 
observations l|SDergel et al.ll2006t ISeliak et alJl2005|) . 

A wide range of dark matter particle candidates have 
been proposed, including axions, the lightest supersym- 
metric particle, and sterile neutrinos. Recently, sterile 
neutrinos have gained renewed interest, since theoretical 
considerations have shown them to be a natural part of 
a minimally extended standard model, the i^MSM. This 
i/MSM, which is the standard model extended with only 
3 sterile neutrinos, ma nages to explain the mass es of the 
active neutrinos ijAsaka. Blanchet fc Shaposhnikovl 
|2005() , the baryon asymmet ry of the uni- 
verse l|Asaka fc Shaposhnikovll2005D. and t he abundance 
of dark matter ijPodelson fc WidrowlllQQ^ . 

Arguably, one of the most appealing aspects of ster- 
ile neutrinos as dark matter candidates is, that they 
will very likely either be detected or rejected dur- 
ing the next couple of years. This is because the 
sterile neutrinos decay virtually at rest into a pho- 
ton and an active neutrino producing a sharp de- 
cay line at E = m s/2 where m.^ is the mass of 
the sterile neutrinos (Polgov fc HansenI l20nl . This 
photon line can now be searched for in cosmic high 
energy data of the background radi ation (for recent 
studie s see iManelli fc Ferraral l|200,'iD : ^oyarsky ct al. 
l)2005)), towards massive structures l i ke galaxy clus- 
ters llAbazaiian. Fuller fc Tucker! 120041: iBovarskv et all 
l200(i() . or towards dark matter structures with a very 
low fraction of baryons (jHansen et alJl2(l?)^ . 

The sterile neutrino was originally introduced in 
order to alleviate discrepancies between the theory 
of cold dark matter structu r e formation and obser- 
vatio ns ijDodelson fc Widrowl Il994t iDolgov fc HansenI 
|2002|) . which indicated that the dark matter particle 
should be warm. Traditionally "warm" means that 
the mass of the dark matter particle should be in 
the keV range, in order to suppress the formation of 
small scale structures. While some of these problems 
have found other solutions, sterile neutrinos have found 
other uses, in cluding an explanation for the pulsar pecu - 
har velocities l)Kusenko fc Segrdl997HFuller et alJ2003j) . 



synth esising early star formation ijBiermann fc Kusenkd 
120061) . and, as mentioned above, as an explana- 
tion for the masses of the active neutrinos and the 
baryon asymmetry |A saka. B lanchet fc Shaposhnikovl 
l2005t lAsaka fc Shaposhnikovll200 5'). Also constraints on 
doube /3 decays have been derived (Bcr uzko v 2005). 

We will here follow the idea of iHansen et al.l l|2002l) 
and look for a decay line towards a region with very little 
baryons. Specifically we will study X-ray "blank" regions 
in the sky with no known X-ray sources. By doing so 
we will search for a signal from the dark matter halo 
of our own Galaxy. We will see that the signal is very 
low for photons with energies of the order 0.1-10 keV, 
which allows us to exclude part of the sterile neutrino 
parameter space. 




Energy (keV) 

Fig. 1. — The blank sky ACIS-S3 spectrum - i.e. a view through 
the Milky Way halo. Also shown is the maximal single gaussian 
emission line from decaying dark matter at a specific energy as 
discussed in Section 



2. X-RAY DATA ANALYSIS 

We have analysed a 0.3-9 keV spectrum extracted from 
many combined Chandra ACIS-S3 blank sky observa- 
tions^ (shown in FIG. [U with CIAO 3.3^ and fitted 

^ http:/ /cxc. harvard.edu/contrib/maxim 
^ http://cxc.harvard.edu/ciao 
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spectral models with the spectral fitting package Xspec 
llArna,iidl ll99(l). 

The hne features at 1.74 keV, 2.1-2.2 keV, and 
7.48 keV are Si ivTa, the Au Maf3 complex, and Ni Ka 
respectively, originating from fluorescence of material in 
the telescope and focal plane (i.e. they are seen in spec- 
tra obtained when ACIS-S3 is stowed and not pointing at 
the sky'^). However, a decay line from dark matter could 
"hide" under these prominent lines. In this case the de- 
cay line should be redshifted by a factor of 1 -I- z where 
z is the redshift of the dark matter emitter. Dark mat- 
ter concentrations located at different distances should 
thus give rise to lines at different observed energies. In 
order to test for this, we extracted spectra from ACIS-S3 
data of two galaxy clusters, A383 and A478, at redshifts, 
z = 0.1883 and 0.0881, respectively. We optimized the 
ratio between expected dark matter signal and emission 
from hot intracluster gas by extracting spectra from the 
cluster outskirts. The X-ray emission from hot intra- 
cluster gas is proportional to the gas density squared, 
whereas the emission of photons from dark matter is di- 
rectly proportional to the dark matter density. Hence, 
in the cluster outskirts the gas surface brightness falls off 
with projected cluster-centric distance, b as, Sgas oc 
(assuming a typical /3-model) whereas the dark matter 
surface brightness is expected to fall off as Sdm oc 
(assuming an isothermal sphere). Specifically, we have 
extracted spectra from an outer radial bin, defined by 
optimizing the ratio of matter in the field of view (pro- 
portional to the signal from decaying dark matter) to 
the X-ray emission from the intracluster gas ("noise") 
described by a /3-profile. The outer radius was chosen as 
the radius with the optimal ratio (which is very close to 
the edge of the ACIS-S3 chip) and the inner radius was 
chosen so the signal to noise ratio had a value of half the 
optimal value. No obvious line features was detected at 
the corresponding redshifted energies'*. 

We then proceeded to firmly constrain the flux from 
any dark matter decay line in the full energy interval 
spanned by the blank sky spectrum. The blank sky spec- 
trum was well fitted (reduced = 1.1 for 540 d.o.f.) 
by a composite model consisting of an exponential plus 
a power law plus four gaussians for the most prominent 
lines. A hypothetical mono-energetic emission line in the 
spectrum was represented by a gaussian, centered at the 
line energy, and with a width, a, given by the instru- 
ment spectral resolution: a k, 0.1 keV for 0.3 keV < E < 
6.0 keV and a « 0.15 keV for 6.0 keV < £; < 8.0 keV^ 
(The velocity of the halo dark matter is producing neg- 
ligible line broadening, w/c ~ 10"'^ — 10~*). For each 
energy in steps of 0.05 keV, we defined a gaussian with 
the instrumental width and maximum at the model value 
of the fit to the broad-band spectrum. The flux of this 
gaussian was then calculated providing an upper limit of 
emission from decaying dark matter at the given energy. 
This is a model independent and very conservative anal- 
ysis, since any X-ray emitters within the field of view, 
and background, contributes to the extracted spectrum. 

^ http:/ /cxc. harvard.edu/proposer/POG/html/nodel. html 
* A fter finisihing this analysis a related paper ap- 
peared JBovarskv et al.l2006l) where an analysis using exactly this 
method of looking at the outer cluster region was made, 
http:/ /acis. mit.edu/acis/spect/spect. html 



If we instead only allowed the decay signal from poten- 
tial decaying dark matter to produce a bump above the 
"base line" spectrum, then we could imp rove the bounds 
by a large factor IBovarskv et al.ll2005|) . However, this 
approach would depend sensitively on the robustness of 
modelling background and emissions from sources con- 
tributing to the "base line" spectrum. 

3. DECAY RATE OF STERILE NEUTRINOS IN 
THE MILKY WAY HALO 

The amount of dark matter from our Galaxy within the 
observed fleld of view is only a minute fraction of the total 
mass of the halo. We use the model for the M ilky Way 
presented in lKlvnin. Zhao fc Somervilld l|2002D . namely 
an NEW dark matter profile with Myir = lO^^M©, virial 
radius ryir = 258 kpc, concentration c = 12, and solar 
distance Rq — 8 kpc. We then integrate out through the 
Milky Way halo to find the expected signal. Changing 
the details of this model does not significantly aflFrect our 
results (e.g. changing the inner density slope between —1 
and zero, or changing the outer slope between —3 and —4, 
changes the predicted signal by less than a factor of 2). 
There is also an uncertainty, possibly as large as a factor 
of 2, arising from using the value of M^r = IO^^Mq. 

The blank sky spectrum was extracted from the full 
angular opening of the Chandra ACIS-S3 chip which is 
8.4' corresponding to a mass within the field of view of 
Mfov — 1-5 X at a mean luminosity distance of 

i^L = 35 kpc, derived from the model mentioned above. 

The luminosity from decaying sterile neutrinos is C = 
E^NT^, where E^ is the energy of the photons, N is the 
number of sterile neutrinos particles in the field of view 
(N = Mfov/rus) and is the decay rate. The flux at a 
luminosity distance, Dl, is; 



F 



C 



4nDl 



EyNT^ 

1^ 



(1) 



The observed flux, F, at a given energy {E^ — ms/2) 
yields a conservative upper limit on the flux from decay- 
ing sterile neutrinos, so eq. Qcan be rewritten as: 
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From the blank sky flux, an upper limit on the de- 
cay rate has been derived as function of line energy and 
plotted in FIG. El 

4. CONFRONTING THE MODEL WITH DATA 

The neutrinos can decay via various channels, and 
the lifetime of the sterile neutrinos is given by 
(jBarger. Phillips &: SarkallT995t iDokov et al]l2000(l : 
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sin (2^) is the mixing angle with the active neutrino, and 
/(ttIs) takes into account the open decay channels so that 
for nis < IMeV, where only the neutrino channel is open, 
/(m,) = 0.86. 

The mixing angle is given by ijDolgovl I2002L 
eqs. 209, 210): 
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Fig. 2. — The decay rate as function of photon energy. The soUd 
Hne is the upper hmit from the Chandra blank sky observations 
viewing out through the Milky Way halo. 

where {Tproduced) is the number of relativistic degrees 
of freedom at the temperature at which the sterile neu- 
trinos are produced. For neutrino masses of « keV, 
the ne utrinos are produced near the QCD phase tran- 
sition l|Dodelson fc Widrowl I1994D. making the distr i- 
bution somewhat non- thermal fPolg ov fc Hansenll2002D . 
The value of g<, for such neutrinos is very conserva- 
tively considered to be between 10.75 and 20, depend- 
ing on the details of the QCD phase transition. We 
will here use g*(Tproduced) = 15 as a reference value. A 
numerical calculation including the details of the QCD 
phase-transition ( Abazaiian 20051) has confirmed that 
the choice — 15 gives good agreement with the an- 
alytical results presented bv ■Dolgov ( 2002^) . S is a free 
parameter taking into account that additional entropy 
may be produced after the sterile neutrinos ha ve been 
created (jAsaka. Kusenko"fc Shaposhnikov 2006() . This 
parameter was sugges ted to be in the range be tween 1 
and 100 ((Asaka. Kusenko fc ShaDOshnikovll200(itl . A is a 
constant depending on which active neutrino the ster- 
ile neutrinos are assumed to mix with. It takes the 
values Ase = 6.7 x 10~® for i^g mixing with z/g and 



As^ = 4.8 X 10- 



for Vfj^.r- The sterile neutrinos are as- 



sumed to account for all dark matter, i.e. floM = 0-30, 
and the Hubble parameter at the present time, Hq, given 
a s h = gn/ (100km/s/ Mpc) 

iBarger. P hillips fc S arkad l|1995f l derived the branch- 
ing ratio for the radiative decay, 1^0+7, to be 
T^/Ttot = (27Q!)/(87r) « 1/128. Now the lifetime can be 
rewritten as: 
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An upper limit on the decay rate is constrained from 
eq. 121 and 

r = ^ « . (6) 



i28r^ 



The lower limit on the lifetime from the data and the 
model predictions are plotted in FIG.|51for various values 
of S and g*. 



We see from FIG. 13 that the blank sky data (solid red 
line) is at least one order of magnitude less restrictive 
than the simplest sterile model predictions (hatched re- 
gion). Models with significant entropy production, S ~ 
30, are excluded in the mass range 4.2 keV < m < 7keV. 
For S > 100 we exclude 2 keV < m < 16keV. 

These constraints are more restrictiv e than the ones 
obtained for the Coma cluster periphery ijBovarskv et alJ 
l200fi|) for masses less than approximately 8 keV. 



= 10-20, S=1 

g=15, S = 5 

0=15, S = 30 

g=15, S=100 




2 4 6 8 

Photon energy [keV] 

Fig. 3. — The lifetime constrained from the flux of the Chandra 
blank sky data (red) and A383 (green). The vMSM prediction for 
5 = 1 and g» = 10 — 20 (hatched) and several variations of S and 
g, (black) have been overplotted. 



5. OTHER CONSTRAINTS 



The simplest 
from below, m, 
llViel et ahl Hcm 
according to a 



models with S* = 1 are bounded 
> 2 keV, using Lya observations 
, and possibly e ven w.^ > 1 4 keV 
recent analysis l|Seliak et al.l l200(i() . 
Th ese constraints are weakned when allo wing for S > 
1 ijAsaka. Kusenko fc Shaposhnikov! [20061) . since addi- 
tional entropy dilutes the momentum of the sterile neu- 
trinos. 

Upper limits on the sterile neutrino mass are de- 
rived fr om the flux of the diffuse photon back- 
ground l)Dolgov fc Hansenl 12002') . and are of the order 
m, < 15 keV for S = 1 (MapcUi fc Ferrara 200^ 
iBovarskv et al.l l2005j) . These bounds are strengthened 
for 5' > 1, since the sterile neutrinos then will decay 
faster. 

A strict lower bound on the sterile neutrino mass 
arises since Big Bang nucleosynthesis only allows ap- 
proximately 0.3 additional relativistic neutrino species 
to be populated at ~ MeV temperatures. This lower 
bound is approximately rus > 50 eV. The Tremaine- 
Gunn bound applied to dwarf speroidals give > 0.5 
keV (Tremaine & Gunn 1979; Lin fc Fabcr 1983). 

The strong claim of an upper m ass limit presented in 
lAbazaiian. Fuller fc TuckeJ l)2004) . under-estimated the 
flux from Virgo cluster by 2 ord ers of magnitude, and is 
unreliable ijBovarskv et alJf2006D . 

Comparing these existing bounds with our findings, we 
conclude that a broad band in neutrino parameter space 
is still open as sho wn in FIG. 01 Taking the results of 
ISeliak et alJ l)2006(l into account, then only masses above 
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approximately 5 ke V are allowed. Ignoring the results of 
iSeliak et alJ pOOfiD . then another band opens in parame- 
ter space, for small masses and large entropy production. 
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Fig. 4. — The rria — S parameter space with constraints from the 
Chandra blank sky data (red, fat, dia gonal, th is paper), Ly-a obser- 
vations (blue, diago nal, conserva t ive [ yiel et al. 12005), and green, 
vertical, ambitious ISeliak et al.l 120061) ). Tremaine-Gunn bound 
m > 0.5 keV (black, horizontal). 



6. CONCLUSIONS 

We have used Chandra blank sky observations, i.e. 
viewing out through the Milky Way halo, to search for a 
possible decay line from dark matter particles. No obvi- 
ous decay line was identified. 

The blank sky 0.3-9 keV X-ray flux is used to constrain 
the parameter space of sterile neutrinos, predicted to de- 
cay at rest to a detectable photon and an active neutrino. 
We find that the entropy production must be limited. 
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S < 100 for the mass range, 2 keV < TOs < 16 keV, and 
even 5" < 10 for masses near = 10 keV. 

This leaves a wide allowed space in parameter space, 
and the sterile neutrino is therefore still a viable dark 
matter candidate. If all Lya constraints are considered, 
then it appears that the sterile neutrino is either rather 
massive, or its momenta are diluted by the entropy pro- 
duction, rendering it a cold dark matter candidate. If 
some of the Lya constraints are relaxed, then there is a 
broad allowed band for smaller masses where the sterile 
neutrino is still a warm dark matter candidate. 

Our analysis is conservative and assumes that the en- 
tire blank sky signal at a given energy is arising from a 
decaying neutrino. By looking for decay lines above the 
expected blank sky emission, the results obtained here 
may be improved significantly. 
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Note added: 

After finishing this paper we became aware of an 
independent analysis, where the possibility of measuring 
decaying dark matter particles from our Galaxy, using 
the blank s ky observation s with XMM-Newton, was 
considered ijBovarskv et all l2006f) . The conclusions 
reached in that paper are qualitatively similar to ours. 

Facilities: Chandra(ACIS) 
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